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ON  THE  DETERMINATION  OF  VERTICAL  SETTLING  OF  THE  BOUNDARY 
OF  A  I-IALF  SPACE  IN  THE  CASE  OF  PLANAR  DEFORMATION 
(A  Dynamic  Problem) 

-y 

N.  M.  Borodachev 


I.  WE  USE  AN  ORTHOGONAL  SYSTEM  OF  COORDINATES  (x,y,z)  THE  y-AXIS.  OF  WHICH 
IS  PERPENDICULAR  TO  THE  BOUNDARY  OF  A  HALF- SPACE. 

Lame's  equations  for  the  elastic  isotropic  half- space  in  the  planar 

case  have,  in  the  absence  of  volumetric  forces ,  the  following  form: 


.  .  <50  ,  , 

O'  +  v-)  —  +  AV2^  —  P  ^7  > 

,,  ,  ,  <30  ,  d-v 

(/.  +  a)  —  +FV^=P— , 


(1.1) 


where  A,  u  are  Lame's  constants  of  elasticity;  u,  u  are  the  projections  of 


the  displacement  vector  on  the  axes  of  the  coordinates  x,  y; 

rj  __  ou  .  dv 
dx  dy  ’ 


p  is  the  density  of  the  material;  and 


Vs 


cP  ,  <P_ 
dx-  ~  dy?  ' 


The  stress  components  are  connected  with  displacement  components  by 


correlations  for  planar  deformations  v.„Leh  appear  as: 


'  T  O  0U 
5*=/J+2’^; 


ay  =  X  0  -f-  2;j- 


du 

dy 


'  dx) 


1 


where 


l 


F  (*.y.  p)  =■ 


Vs- 


sH 


1-p'+i'*0(X,  y,t)dt. 


The  solution  of  equation  (1.8)  which  is  bounded  at  infinity  appears  as: 
F(a,y,p)—  A  exp  ^[)+  BexP  (_>'|//'“5  +  t)  ’  f1*9) 


where  A  and  B  are  determined  from. conditions  on  the  boundary  of  the  half¬ 
space.  : 

Successively  applying  to  (1.3),  (1.6)  and  (1.7)  a  Laplace  transform 
with  reference  to  the  variable  t,  and  a  Fourier  transform  with  reference  to 
the  variable  x,  we  shall  have  the ' following  relationship  between  the 
Laplace -Fourier  transformants  of  the  functions  v  (x,  y,  t),  ay  ( x,  y,  t) , 
tvy(x,  y,  t),  and  of  the  function  F  (a,  y,  p) .  Substituting  into  these 
relationships  functions  F  (a,  y,  p)  from  (1.9),  and  also  assuming  that 
y  =  0,  we  obtain  the  following  expression  for  the  Laplace -Fourier  transfor¬ 
mant  of  the  functions  v  (x,  0,  t),  cry(x,  0,  t)  and  Txy  (x,  0,  t) : 


(1,10) 

(1.11) 

(1.12) 


II.  OUR  PROBLEM  CONSISTS  IN  DETEFli-HKINC-  THE  FUNCTIONS  v  (x,  0,  t)  IF 
O'.j  (x,  0,  t)  AND  T,y(x,  0,  t)  ARE  KGS. 

VJe  obtain  the  solution  of  this  problem  by  utilizing  the  results  of 

paragraph  I  of  this  article.  In  fact,  when  the  functions  ay (x,  0,  t)  and 

,  (x,  0,  t)  are  known,  then  the  Laplace-Fourier  Transformants  of  the 
yC  j 

fuueuioas  Zy  (&,.  0,  p)  and  ~£  (a,  0,  p) ,  generally  speaking,  also  may  be 


considered  as  known.  On  the  "basis  of  this,  from  (l.ll)  and  (1.12)  egressions 
for  A  (a,  p)  and  B  (a,  p)  can  easily  "be  found.  Substituting  the  expressions 
for  A  (a,  p)  and  B  (a,  p)  into  (1.10)  and  successively  applying  the  inversion 
theorem  for  the  Laplace  and  Fourier  transform  /2/,  we  can-  obtain  the  formula  for 
the  determination  of  the  function  v  (x,  0,  t). 

i  ,  .  .  ' 

Let  us  examine  in  more  detail  the  case  ‘L^y(x,  0,  t)  =  0.  Substituting 
the  expressions  A  (a,  p)  and  B'(a,  p)  from  (l.ll)  and  (1.12)  for  the  case 
(x,  0,  t)  =  0  in  (1.10)  after  algebraic  transformations,  we -obtain  the 
following  formula  for  the  Laplace-Fourier ■ transforms  of  the  function 
V  (x,  0,  t): 


•  ^  (a.  0 ,p)  = 


i,.  (a,  0 ,p) 

”*  6  ~  ~ 
^  II  (p-iac^h ) 

k~\ 


+ (?+§)’ 


where  x^'  are  the  roots .  of  the  equation  ■ 

/  c2  \  7 

xc_8-r1+[24- 16^jx2-|-^16-^ 

The  roots  of  equation  (2.3)  depend  on  the  relation  ’ 

c"i  '  ^  1  —  2v 

Cj  X  -j-  2  —  2v 

where  y  is -the  Poisson’s  ratio. 

Successively  applying  the  theorems  of  inversion  to  (2.l)  for  the  Laplace 
and  Fourier  transforms,  we  obtain: 


P  :)**•/ 


ij+/~ 

Jpf-ioJC 


^v(q,  0,/>)  V 

6 

n  (/>— ««**«) 

k~\ 


-r  (x,  0,  0  = 


The  final  stage  of  solving  the  present  problem  is  the  calculation  of 
the  contour  integrals  entering  into  expression  (2.3).  However,  for  practical 
use  it  is  better  to  have  a  formula  for  function  v  (x,  0,  t)  into  which  real 
integrals  -  not  contour  integrals  -  enter.  For  obtaining  such  a  formula, 
we  apply  successively  to  (2.1)-)  the  theorems  about  the  convolutions  for  the 
Laplace  transform  in  reference  to  the  variable  t,  and  for  the  Fourier  trans¬ 
form  in  reference  to  the  variable  x  /2 /.■  As  a  result  we  get: 


where 


t 

■  v(x,0,t)=  —  [d-  f  V  .  (2.4) 

.  J  J  01  • 


«•(«,*) 


:  0  at 

'  c'  ■  ; 

.  A«f*  k)+/4s(?2— . 

**  cT1  . 

Ao  ( fi  -rf^-rA,  (?1 

■at  ’  c,  ^ ‘  <  x3c,  ’ '  . 

Ao  ( ft  ~r.f 2)  4"  A]  (?i  ~.bj)  +  A-,  (c2  —  ti,)  -j-  A3(oj  -j-  <!)j) 

•  ■  •  . . ,  ‘,<jl 

an.  ^ 


(2.5)' 


H-f -1 


7.3C2 


(-if 


/•=1±_L.  j  i. 

•  4  ■  (f-f) 


A,  : 


S  <3  , 

■  H 


f  H-4) 


•;  a,= 


(-if.  M)V-H? 


(2.6) 


H-4) 


/>  (i  *)  =  In 
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(2.7) 
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Thus,  a  general  solution  to  the  problem  of  determining  the  vertical 
components  of  the  displacement  vector  of  the  boundary  of  the  half-space  in 
a  planar  case  is  reached  by  formula  (2.4)  when  t  >  0  stress  acts  on  the 
boundary  of,  the  half-space  at  Oy  (x,  0,  t.)and 't ^  (x,  0,  t)  =  0. 

The  obtained  solution  is  valid  when  the  following  conditions  are  fulfilled 

a.  At  t  ^0  the  half-space  is  in  the  state  of  rest; 

b.  r/y  (x,  0,  0)  =  0; 

c.  All  roots  of  the  equation  (2.2)  are  simple  and  real; 

d.  When  the  inequalities  are  valid  0  <  c2  <  c2<c,  <xsC2  <%i  c2> 

where  •/,,  yj,  y.s  are  the  positive  roots  of  the  equation  (2.2).  For  instance, 
at  y  -  0.25  conditions  8  and  2  are  fulfilled.  However,  in,  case  of  necessity 
this  method  will  give  a  corresponding  solution,  also,  when  the  above  enumerated 
conditions  are  valid. 

III.  IN  THE  FOLLOWING  IS  SHOWN  AN  APPLICATION  OF  THE  GENERAL  FORMULA  (2.4) 

FOR  THE  SOLUTION  OF  INDIVIDUAL  PROBLEMS. 
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av  (X,  0,  f.)=~Po(x)H(t), 


Problem  a,  Let 


(3.  1) 


where  o  (:•:)  is  the  Dirac  delta- function,  and  H  (t)  is  a  discontinuous 
function  determined  by  the  correlation 


Ht=.  (0  ctKO, 
U  at^>0. 


Phe  case  under  examination  is  graphically  depicted  in  Pig.  1. 


In  this  case 

AcJx.OJ)  =  —  P 3  (x)  o  (t) 
dt  ■ 

and  formula  (2.4)  gives 


„  (X,  0,  t)  =  j  dr.  j  v  (*,  t)Ht-  W  (*  -  5)  dl. 

0  — 3 


Utilising  the  properties  of  the  delta- function  /k/  we  finally  get 


«lje,0,  t)  =  —V(x,  t), 


(3.2) 


where  the  function  0  is  determined  by  the  formulas  (2.5)- (2.8),  wherein  the 
variables  £  and  't  are  substituted  by  x  and  t  respectively..  When  the  boundary 
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value  is  calculated  "by  the  function  v  (x,  0,  t)  determined  by  the  expression 
(3-2)  at  t  — ■>  03 ,  then  we  arrive  at  a  known  formula 


lim  v(x,  0,  t) 

t .. 


nE  U| 


(3.  3) 


which  is  obtained  from  the  investigations  of  Flamant  /5/,  and  which  corres¬ 
ponds  to  the  case  of  a  statical  application  of  stresses  Gy  (x,  0)  =  -P  6  (x). 
Here  the  expression  E  =  2  (l  +  v)  /i  is  the  modulus  of  the  normal  (perpendi¬ 
cular)  stress.  It  must  be  noted  that  formula  (3.3)  does  not  determine  the 
absolute  value  of  the  vertical  component  of  the  displacement  vector  of  any 
arbitrary  point  located  on  the  boundary  of  the  half- space,  but  the  differ¬ 
ence  of  the  vertical  components  of  the  displacement  vector  at  the  point  with 
the  abscissa  (x)  and  a  certain  point  also  located  on  the  boundary  but  removed 
at  a  distance  of  (d)  from  the  origin  of  the  coordinates. 

Problem  b a  let 

=  (3.4) 

The.  examined  case  is  graphically  depicted  in  Fig.  2.  In  the  given  case 


—  0  (x,  0,  t)  —  —  q  H  (d—  \x\)  S  (t) 
dt  v 

and  formula  (2.5)  gives 

t  x+a 

v(x,0,t)= 

0  x—a 

Changing  the  order  of  integration  in  this  expression  and  accomplishing  the 
integration  by  'C  ,  we  get 

v  (x,  0,  0  =  “  J  W  (£,  t)  & .  (3.  5) 

x—a 
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Integration  in  formula  (3.5)  can  be  accomplished  without  undue  difficulties. 
But  since  the  final  expression  in  this  case  for  function  v  (x,  0,  t)  is  too 
cumbersome ,  we  will  not  write  it  out. 
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